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SPONTANEOUS  REACTIVATION  OF  ACETYLCHOLINESTERASE  FOLLOWING 
ORGANOPHOSPHATE  INHIBITION.  III.  STUDIES  WITH 
/J-NITROPHENYL  METHYLPHENYLPHOSPHINATE 


1.  INTRODUCTION. 

A large  number  of  papers  in  recent  years  have  been  devoted  to  the  problems  of  the 
biochemical  mechanism,  toxic  action,  and  therapy  involving  organophosphorus  compounds  which 
are  inhibitors  of  cholinesterases.  Regarding  therapy,  the  development  of  oximes  as  antidotes 
against  cholinesterase  poisoning  followed  the  observation  that  the  very  slow  restoration  of  enzymatic 
activity  obtained  with  certain  phosphorylated  cholinesterases  was  accelerated  in  the  presence  of 
nucleophilic  agents  such  as  choline  and  hydroxylamine.1 ,2  The  use  of  the  latter  nucleophiles  was 
followed  by  hydroxamic  acids  and  additional  ketoximes  and  aldoximes.  However,  even  the  present- 
day  oximes  such  as  2-PAM*  and  TMB-4**  are  of  limited  value  in  the  therapy  of  organophosphate 
poisoning.3,4 

Our  interest  in  the  spontanebus  reactivation  of  phosphorylated^  acetylcholinesterase 
( AChE,  EC  3.1 .1 .7)  is  directed  at  understanding  the  basic  parameters  and  mechanism  of  this 
phenomenon.  Such  knowledge  is  mandatory  if  one  hopes  to  advantageously  perturb  this  system 
as  part  of  a therapeutic  approach  in  organophosphate  poisoning. 

Our  first  studies  in  this  area  were  carried  out  using  a series  of  paru-substituted  phenyl 
methylphosphonochloridates.5  These  spontaneous  reactivation  studies  produced  the  first  evidence 
for  hydrophobic  interactions  participating  in  the  spontaneous  dephosphonylation  of  inhibited  AChE. 
We  have  subsequently  reported  the  first  kinetic  studies6  on  the  spontaneous  reactivation  of  AChE 
inhibited  with  Sarin. '' 

In  support  of  our  continuing  effort  in  the  area  of  spontaneous  reactivation  studies,  we 
synthesized  p-nitrophenyl  methylphenylphosphinate  (NMPP).  The  reasons  for  the  selection  of  a 
phosphinate  in  our  current  studies  were  twofold.  First,  in  the  case  of  phosphinate-inhibited  AChE, 
aging  should  be  excluded  from  reaction  considerations  since  a susceptible  carbon-oxygen-phosphorus 
bond,  necessary  for  aging,  is  absent.  Second,  an  additional  positive  aspect  resulting  from 
phosphinate  use  concerns  the  amino  acid  functional  in  acid-catalyzed  aging.  Since  there  is  evidence 
that  the  same  amino  acid  is  involved  in  substrate  conversion,7 ,8  the  use  of  phosphinates  precludes 
an  approach  where  the  action  of  this  amino  acid  may  be  affected.  Consequently,  the  spontaneous 
reactivation  process  is  effectively  isolated  and,  therefore,  approaches  to  its  advantageous  perturbation 
may  be  investigated. 

Our  work  to  date  with  NMPP  supports  the  rationale  that  phosphinates  are  ideally  suited 
to  spontaneous  reactivation  studies.  Of  special  interest  is  the  fact  that  this  report  also  documents 
the  first  spontaneous  reactivation  studies  of  AChE  following  inhibition  by  a phosphinate. 


* Pyridinium  aldoxime  methochloride. 

**  1 ,l'  Trimethylcne-A(j(4-formylpyridinium  bromide)dioxime. 

t The  term  phosphorylation  will  be  used  to  include  phosphorylation,  phosphonylation,  and  phosphinylation. 
[Steinberg,  G.  M.,  Lieske,  C.  N.,  Boldt,  R.,  Goan,  J.  C.,  and  Podall,  H.  E.  Model  Studies  for  the  Reactivation  of 
Aged  Phosphonylated  Acetylcholinesterase.  Use  of  Alkylating  Agents  Containing  Nucleophilic  Groups.  J.  Med. 
Chem.  13, 435  (1970)) . 
tt  Isopropyl  methylphosphonofluoridate. 
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11.  MATERIALS  AND  METHODS. 


A.  General. 

The  zwittcrionic  3-(N-morpholino)propanesulfonic  acid  (MOPS)  and  N.N-fc/M 2-hydroxy- 
ethyl  Jglycine  (Bicine)  buffers  were  described  by  Good  el  al?  and  are  available  commercially  from 
several  sources.  MOPS  was  used  in  the  pH  7.60  studies;  and  Bicine,  in  the  pH  9.10  studies.  Each 
was  0.10  M and  contained  0.01  M Mg+  , 0.002%  NaNj,  and  0.01%  bovine  serum  albumin. 

The  0.134  M phosphate  buffer  of  pH  6.94  was  prepared  from  0.134  M K2HPO4  and 
0. 1 34  M KH-.PO4  solutions.  Sodium  azide.  20  mg/1,  was  added  to  prevent  bacterial  growth.  A 
1 1 (v/v)  dilution  of  this  buffer  with  double  distilled  water  gave  a pH  of  6.90  for  the  resulting 

O. 067  M phosphate  buffer.  This  buffer  was  used  in  the  pH  6.90  studies.  All  pH  values  were 
determined  at  25.0°C  using  a Beckman  research  pH  meter,  Model  101900,  equipped  with  a Markson 
No.  1808  Polymark  combination  electrode. 

The  eel  acetylcholinesterase  (AChE,  EC  3. 1.1. 7)  was  purchased  from  Worthington 
Biochemical  Corporation.  It  was  obtained  as  a purified,  salt-free  powder,  4000  units/mg.  Two  and 
one-half  milligrams  of  this  powder  was  dissolved  in  1 .5  ml  of  a previously  boiled  solution  (pH 
approximately  7.4)  containing  0.225  M KC1,  0.10%  gelatin,  and  0.020%  sodium  azide  to  give  the 
enzyme  concentrate.  The  addition  of  5.0  /al  of  a fiftyfold  dilution  of  the  concentrate  to  2.00  ml  of 
010  M MOPS  buffer,  pH  7.60,  containing  5 0 /al  of  1 .58  M phenyl  acetate  in  acetonitrile,  gives  an 
absorbance  change  at  272-5  nm  of  ca.  0.60  absorbance  unit/min  (25.0°C).  Assay  concentration  of 
the  substrate  equals  4 X 10“^  M. 

B.  Preparation  of  ;>-Nitrophenyl  Methylphenylphosphinate. 

Methylphenylphosphinyl  chloride  (550  mg,  3.15  mmoles)  dissolved  in  5 ml  of  benzene 
was  added  in  drops  over  a 10-minute  period  to  a magnetically  stirred  solution  of  p-nitrophenol 
(414  mg,  2.98  mmoles)  and  triethylamine  (290  mg,  2.87  mmoles)  in  10  ml  of  benzene.  Upon 
completion  of  the  addition,  the  stirred  reaction  mixture  was  maintained  at  room  temperature  for 
several  hours.  It  was  transferred  to  a separatory  funnel  and  washed  progressively  with  small  aliquots 
of  pH  4 acetate  buffer  (3X),  pH  7.5  MOPS  buffer  (3X),  and  finally  with  distilled  water  (2X).  After 
the  nonaqueous  layer  was  dried  over  MgSO^,  the  benzene  was  removed  in  vacuo,  without  heating, 
to  give  a viscous  yellow-orange  oil.  The  oil  solidified  on  standing  overnight.  Two  recrystallizations 
trom^ethyl  ether  gave  1 58  mg  (19%  yield)  of  product.  The  very  light  tan  crystals  melted  at  85°  to 
86.5°C  (uncorrected).  The  integrated  NMR  spectrum  was  in  consonance  with  the  assigned  structure 
and  showed  no  evidence  of  impurities.  Anal.  Calc  for  C13H.2O4NP:  C,  56.33;  H,  4.36;  N,  5 05 

P,  1 1 . 1 7.  Found:  C,  56.4;  H,  4.4;  N,  5.0;  P,  1 1 .2. 


The  mass  spectrum  of  NMPP  was  obtained  on  a Hitachi  Model  RMU6E,  single-focusing, 
magnetic  deflection  mass  spectrometer.  The  solid  sample  was  placed  near  the  source  of  the 
spectrometer,  using  the  solid  probe  feature  of  the  instrument.  The  tip  of  the  probe  was  surrounded 
by  an  electrically  heated  copper  block;  its  temperature  controlled  the  vaporization  of  the  sample. 

1 he  ionizing  voltage  was  70  volts.  A list  of  the  relative  intensities  obtained  with  the  sample  heater 
held  at  200°C  appears  on  page  13. 


D. 


Hydrolysis  of  p-Nitrophenyl  Methylphenylphosphinate. 


1 


The  hydrolysis  of  NMPP  in  phosphate  buffer  was  monitored  spectrophotometrically  at 
400  nm.  The  reaction  was  initiated  by  introducing,  with  rapid  mixing,  6.0  pi  of  a 4.87  X 10'^  M 
solution  of  NMPP  in  acetonitrile  into  an  UV  cell  containing  3. 00  ml  of  0.067  M phosphate  buffer 
of  pH  6.90  at  25.0°C.  The  data  obtained  were  calculated  using  equation  1 . Linear  first-order 
kinetics  were  observed  for  the  1 .9  half  lives  monitored.  At  completion  of  the  hydrolysis,  the 
production  ofp-nitrophenol  was  stoichiometric.  The  results  of  three  runs  under  the  cited  conditions 
gave  a k^j  of  0.0219  min'*  ± 2.59%,  corresponding  to  a half  life  of  31 .64  minutes. 

The  hydrolysis  of  NMPP  in  MOPS  buffer  was  monitored  spectrophotometrically  at 
400  nm.  The  reaction  was  initiated  by  introducing,  with  rapid  mixing,  4.0  pi  of  a 5.27  X 10'"  M 
solution  of  NMPP  in  acetonitrile  into  an  UV  cell  containing  3 00  ml  of  0.10  M MOPS  buffer  of 
pH  7.60  at  25.0°C.  The  data  obtained  were  calculated  using  equation  1 . Linear  first-order  kinetics 
were  observed  through  the  1.97  halflives  monitored.  At  completion  of  the  hydrolysis,  the 
production  of  p-nitrophenol  was  stoichiometric.  The  results  of  three  runs  under  the  cited  conditions 
gave  a k0jjS£j  of  0.0066  min'*  ± 1.8%,  corresponding  to  a half  life  of  105.0  minutes. 

L.  Inhibition  of  Lei  Acetylcholinesterase  by  p-Nitrophenyl  Methylphenylphosphinate. 

The  inhibition  of  AChE  by  NMPP  was  studied  by  an  adaption  and  modification  of  the 
Hart  and  O’Brien1  0 method  for  obtaining  the  inhibition  constants.  The  kinetic  experiments  were 
initiated  by  mixing  reactants  by  an  Aminco-Morrow  stopped-flow  apparatus  which  was  attached  to 
a Beckman  DUR  monochromator.  A data  acquisition,  storage,  and  retrieval  system  (DASAR, 
American  Instrument  Company)  was  used  for  visual  examination  and  acquisition  of  the  data.  A 
high-speed  paper-tape  punch,  interfaced  with  the  DASAR,  was  employed  to  produce  a permanent 
record  of  experimental  results  on  paper  tape.  The  permanent  data  record  on  paper  tape  was 
transferred  to  the  Univac  1 108  computer  files.  Program  and  data  files,  as  described  by  Horton  etal.,1 1 
were  called  and  manipulated  by  use  of  a remote  terminal  keyboard. 

The  reactant  solutions  were  prepared  at  twice  the  desired  final  concentration  after  mixing 
in  the  stopped-flow  instrument.  The  AChE  solution  in  0. 134  M phosphate  buffer  was  prepared  by 
diluting  10.5  pi  of  enzyme  concentrate  (section  II,  A)  with  30  ml  of  the  buffer.  Final  buffer 
concentration  was  0.067  M after  mixing  and  the  pH  was  6.90.  Just  prior  to  use,  solutions  of  (1 ) the 
substrate,  p-nitrophenyl  acetate,  and  (2)  the  inhibitor,  NMPP,  plus  substrate  were  prepared  in  4% 
ethanol/96%  distilled  water  (v/v).  This,  of  course,  resulted  in  a 2%  ethanol  solution  in  each  case 
following  mixing  in  the  stopped-flow  apparatus. 

The  p-nitrophenyl  acetate  solution  was  used  to  determine  control  velocities,  and  the 
p-nitrophenyl  acetate  plus  NMPP  solution  was  used  for  the  inhibition  runs.  The  concentrations  of 
substrate  were  identical  in  the  control  velocity  runs  and  the  inhibition  runs.  In  all  determinations, 
the  final  substrate  concentration  was  1.00  X 10  3 M;  and,  in  the  inhibition  determinations,  the 
final  concentration  of  NMPP  was  2.67  X 10*6  M.  The  temperature  of  the  reaction  cell  in  the 
stopped-flow  apparatus  was  controlled  at  25°C  by  circulating  water  from  a constant  temperature 
bath.  Residual  enzymatic  activity  and  background  hydrolysis  of  the  substrate  after  complete 
inhibition  were  determined  and  were  incorporated  in  the  calculations  of  the  inhibition  constants 
resulting  from  five  determinations  (see  page  14). 
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F.  Spontaneous  Reactivation  of  Methylphenylphosphinylated  Eel  Acetylcholinesterase. 

A detailed  description  of  the  technique  involved  in  a spontaneous  reactivation  study  of 
methylphenylphosphinylated  AChE  at  pH  7.60  is  given  below.  Other  kinetic  investigations  at 
pH  7.60,  as  well  as  those  at  pH  9.10,  were  handled  similarly. 

An  enzyme  stock  solution  was  prepared  by  adding  10.0  n\  of  enzyme  concentrate 
(section  II,  A)  to  1.5  ml  of  0.10  M MOPS  buffer  of  pH  7.60.  For  a control,  0.500  ml  of  this  stock 
solution  was  added  to  49.50  ml  of  the  buffer.  Phenyl  acetate  assays  of  this  1 to  100  dilution  showed 
an  activity  of  0.566  absorbance  unit/min  at  272.5  nm,  allowance  being  made  for  the  nonenzymatic 
hydrolysis  of  phenyl  acetate  in  the  buffer  medium  (0.002  absorbance  unit/min). 

The  remaining  1 .0  ml  of  enzyme  stock  solution  was  inhibited  with  10.0  #d  of  4. 87  X 10 
M NMPP  dissolved  in  acetonitrile.  Three  minutes  later,  0.500  ml  of  the  inhibited  enzyme  solution 
was  added  via  a Chromatronix  calibrated  sample  loop  to  a 1 /2-inch  by  6-3/4-inch  column  bed  of 
L'ltrogel  AcA44  (4%  acrylamide  and  4%  agarose)  and  eluted  with  the  buffer  at  a flow  rate  of 
0.50  ml/min.  The  zero  time  for  reactivation  was  taken  as  the  time  when  the  0.500-ml  sample  of 
inhibited  AChE  had  passed  into  the  gel  bed  of  the  column.  The  first  17.00  ml  was  collected  and 
diluted  to  50.0  ml  with  buffer.  During  the  course  of  this  elution,  another  aliquot  of  the  inhibited 
enzyme  was  checked  for  residual  activity  with  none  being  detected.  Previous  tests  established  that 
all  the  AChE  and  no  inhibitor  was  present  in  the  17.00-ml  fraction  collected  from  the  column. 

The  50.0-ml  fraction  was  incubated  at  25.0°C.  At  appropriate  intervals,  1 00-ml  aliquots 
of  this  solution  were  assayed.  Activity  was  calculated  from  the  phenyl  acetate  assays,  the  reaction 
being  monitored  in  this  manner  for  ca.  100  hours.  Corrections  were  made  for  the  loss  of  activity 
observed  in  the  control  during  this  time  period.  The  results  of  this  spontaneous  reactivation 
experiment  gave  k0^s  = 0.02120  hr"^  with  a 100.1%  return  of  activity  (Figure  1).  A summary  of 
our  spontaneous  reactivation  experiments  is  presented  on  page  15. 

G.  Calculation  of  the  Observed  Spontaneous  Reactivation  Rate  of  Methylphenyl- 
phosphinylated Eel  Acetylcholinesterase. 

A modification  of  the  computer  program  originally  developed  by  Hayo  and  Wilcoxson1  2 
was  used  in  these  calculations.  The  FORTRAN  program  (EXPFIT)  Fits  the  reactivation  data  to  the 
exponential  rate  equation 

V = V^fr-e^obs^o)]  0> 

where  V = velocity  function  (e.g.,  zero-order  rates,  absorbances,  and  % reactivation),  value  of 
velocity  function  at  infinite  time  (limiting  value),  k0^s  = observed  First-order  rate  constant,  and 
tg  = the  time  value  when  the  velocity  function  is  equal  to  zero  (abscissa  intercept  of  graphical 
treatment).  An  iterative  least-squares  procedure  is  used  to  Fit  (V,  t)  data,  and  the  best  values  for 
V*.  k0bs,  and  tg  are  calculated.  The  validity  of  equation  1 has  been  previously  reported.5  Calculation 
of  first-order  kinetic  constants  when  the  experimental  end  point  is  not  available,  using  the  Hayo  and 
Wilcoxson  method,  offers  a distinct  advantage  over  the  more  common  procedure  of  Guggenheim.1  3 
The  former  method  does  not  require  that  the  kinetic  parameters  be  spaced  at  constant  time 
intervals.  A complete  documentation  of  program  EXPFIT  is  presented  in  appendix  A. 
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RUN- 1 0 : EXP2854 


TIME  IN  HOURS 


Figure  1 . Spontaneous  Reactivation  of  Methylphenylphosphinylated  Eel  Acetylcholinesterase 

(pH  7.60,  0.10  M MOPS  buffer,  25.0°C)  / 

In  addition  to  calculating  the  kinetic  parameters,  program  EXPFIT  also  creates  a separate 
plotting  parameter  file  to  provide  the  capability  of  examining  the  results  graphically.  This  plotting 
file  can  be  subsequently  read  by  program  EXPFITPLOT  which  provides  the  required  machine 
control  commands  to  operate  the  off-line  CALCOMP  plotter.  Appendix  B contains  a complete 
documentation  of  program  EXPFITPLOT. 

f 

H.  Induced  Reactivation  of  Methylphenylphosphinylated  Eel  Acetylcholinesterase. 

A detailed  description  of  the  technique  involved  in  our  induped  reactivation  studies  of 
NMPP-inhibited  AChE  is  given  below.  All  experiments  were  carried  out  at  pH  7.60  in  0.10  M MOPS 
buffer  (25.0°C).  The  oximes  used  were  2-PAM,  TMB-4,  and  MINA.*/  Each  oxime  study  was  treated 
in  the  same  manner.  j 

I 

An  enzyme  stock  solution  was  prepared  by  adding  lO.Oifi  of  enzyme  concentrate 
(section  II,  A)  to  1 .5  ml  of  0.10  M MOPS  buffer.  A control  was  j/repared  by  adding  0.500  ml  of 
this  stock  solution  to  49.50  ml  of  the  buffer.  Phenyl  acetate  assays  of  this  1 to  100  dilution  showed 


* Monoisonitrosoucetone. 
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an  activity  of  0.577  absorbance  unit/min  (272.5  nm).  The  nonenzymatic  hydrolysis  of  phenyl 
acetate  in  the  buffer  (0.002  absorbance  unit/min)  was  corrected  in  each  assay. 


To  the  remaining  1 .0  ml  of  enzyme  stock  solution  was  added  10.0  Ml  of  4.87  X 10"^  M 
NMPP.  The  inhibition  was  allowed  to  proceed  at  25.0°C  for  3 minutes.  Exactly  0.500  ml  of  the 
inhibited  enzyme  solution  (complete  inhibition  confirmed  in  an  independent  assay)  was  added  to  a 
1 '2-inch  by  6-3/4-inch  column  bed  of  Ultrogel  AcA44  and  eluted  with  the  buffer  at  a flow  rate  of 
0.50  ml/min.  The  first  6.00  ml  of  eluent  were  discarded  and  the  7.00-  through  16.90-ml  fraction 
was  collected  in  a 10.0-ml  volumetric  flask. 

At  zero  time,  100.0  m1  of  lO"*-  M oxime  was  added  to  the  gel-filtered  enzyme  aliquot  to 
yield  a final  volume  of  10.0  ml  (10'^  M in  oxime).  The  10.0-ml  aliquot  was  maintained  at  25.0°C 
and  assayed  in  l hour.  In  the  case  of  MINA,  additional  assays  were  carried  out  during  a 4-hour 
period.  For  all  assays,  95.0  m1  of  the  inhibited  enzyme/oxime  solution  was  added  to  1 .90  ml  of 
0. 10  M MOPS  buffer.  Thus,  oxime  concentration  in  the  assay  solution  was  5.0  X 10'6  M.  Activity  was 
calculated  from  phenyl  acetate  assays.  For  correction  purposes,  an  oxime  control  was  determined 
tor  each  experiment  to  ascertain  the  amount  of  reversible  inhibition  exhibited  by  the  oxime  at 
5.0  X 10  M in  the  assay  aliquot.  In  addition,  oxime-catalyzed  hydrolysis  of  the  phenyl  acetate  in 
the  butter  was  measured  and  the  resulting  value  was  incorporated  in  the  final  calculations.  The 
results  ot  these  oxime-induced  reactivations  are  summarized  on  page  16. 

HI.  RESULTS  AND  DISCUSSION. 

The  NMPP  used  in  these  studies  was  prepared  from  methylphenylphosphinyl  chloride  and 
p-nitrophenol.  The  NMR  of  the  starting  chloride  indicated  purity  >98%.  However,  since  the  mass 
spectrum  ot  the  chloride  indicated  small  quantities  of  a variety  of  compounds  (figure  2),  we 
characterized  the  NMPP  obtained  by  melting  point,  elemental  analysis,  NMR,  and  mass  spectroscopy 
(section  II,  B and  C). 

The  mass  spectrum  of  NMPP  is  given  in  table  1 . The  relative  intensities  were  found  to  be 
dependent  on  the  conditions  for  heating  the  sample,  especially  the  temperature  of  the  probe  at  the 
time  the  sampk  was  introduced.  The  intensities  given  in  table  1 were  obtained  with  the  sample 
heater  held  at  200°C’,  a temperature  which  induced  the  solid  to  volatilize  quickly.  Table  1 also  lists 
the  observed  metastable  peaks.  Metastable  peaks  are  very  broad  signals  which  can  appear  at 
nonintegral  mass  values,  whose  position  is  determined  according  to  the  formula  m = (mi)  ^/mi , 
where  mj  is  the  parent  ion  and  m2  the  daughter  ion  of  the  fragmentation  involved.  The  origin  of 
all  the  strong  peaks  in  the  spectrum  can  be  rationalized  by  cleavages  in  which  no  complicated 
rearrangements  are  required.  In  this  regard,  it  is  interesting  to  note  that  there  is  usually  considerable 
rearrangement  of  hydrogen  atoms  in  the  fragmentation  of  the  organic  phosphates  and  phosphonates 
which  is  absent  in  this  phosphinate. 

The  impurity  peak  at  m/e  294  may  be  the  pyro  compound 

0 O 

II  II 

0— P— O— P— 0 

1 I 

ch3  ch3 

which  was  one  of  the  impurities  detected  in  the  mass  spectrum  of  the  methylphenylphosphinyl 
chloride  used  to  prepare  NMPP  (figure  2).  The  fact  that  all  strong  peaks  can  be  accounted 
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for  without  involving  the  presence  of  impurities,  coupled  with  the  stoichiometry  of  p-nitrophenol 
production  in  the  hydrolysis  studies,  permits  a purity  >99.5%  to  be  assigned  to  the  NMPP  used  in 
these  studies. 


r 


Figure  2.  Structural  Formulas  of  Trace  Impurities  in  Methylphenylphosphinyl  Chloride 


Table  1.  Mass  Spectrum  of  p-Nltrophenyl  Methylphenylphosphlnate 


To  additionally  characterize  NMPP  and,  as  a prelude  to  our  enzymatic  studies  with  this 
compound,  its  hydrolysis  was  examined  spectrophotometrically  at  pH  6.90  in  0.067  M phosphate 
buffer  and  at  pH  7.60  in  0.10  M MOPS  buffer.  The  half  life  in  the  MOPS  buffer  (pH  7.60)  of 
105.0  minutes  compared  to  the  half  life  in  the  phosphate  buffer  (pH  6.90)  of  31 .64  minutes 
demonstrates  the  advantage  of  using  a Good  9 buffer  when  it  is  desired  to  keep  buffer  participation 
to  a minimum:  that  is,  in  spite  of  the  higher  pH  in  which  the  MOPS  hydrolysis  runs  were  carried 
out,  the  rate  is  ca.  3-1/2  times  slower  than  observed  in  the  phosphate  buffer  of  lower  pH.  We  have 
reported  similar  results  with  Bicine9  buffer  in  the  hydrolysis  of  O-methyl  -nitrophenyl  phenacyl 
methylphosphonate  oxime.1  4 

For  additional  comparisons,  one  may  approximate  a kQ^j  from  the  MOPS  hydrolysis 
runs  at  pH  7.60  by  neglecting  the  small  contribution  of  MOPS  to  the  reaction  rate.  Letting 
0.0066  min  * = kQ^l3.98  X 10'^],  kQ^  = ca.  1658  M'*  min'*.  Clearly  NMPP  is  more 
susceptible  to  alkaline  hydrolysis  than  p-nitrophenyl  diethylphosphinate  (kQ^  = 128  M'*  min'*), 
p-nitrophenyl  J/-n-propylphosphinate  (kQj^  = 93.4  M'*  min'*),  or /7-nitrophenyl  <//-«-butyl- 
phosphinate  (kQ^  = 81.5  M'*  min'*).15 

The  inhibition  of  AChL  by  NMPP  is  accompanied  by  the  formation  of  a covalent 
phosphorus-enzyme  linkage  at  the  serine  active  site.  This  inhibition  is  a rather  facile  reaction* 
as  evidenced  by  its  rate  constant,  kj,  ot  2.95  X 104  M'*  sec'*  (table  2).  A similar  value, 
1.32X10  M * sec  * , was  reported  by  Fukuto  ct  al.,  using  fly-brain  cholinesterase  and  p-nitrophenyl 
ethylphenylphosphinate.1  6 


Table  2.  Constants  for  the  Inhibition  of  Eel  AChE  by  NMPP.  Substrate, 
p-nitrophenyl  acetate;  pH  6.90,  0.067  M phosphate  buffer;  25°C 


Run  No. 

Kd  X 10* 

k2 

X 10‘4 

M 

-1 

sec 

-1  -1 
M sec 

1 

7.74 

0.230 

2.97 

2 

8.29 

0.258 

3.11 

3 

6.22 

0.175 

2.82 

4 

8.74 

0.268 

3.06 

5 

7.15 

0.198 

2.78 

Mean 

7.63 

0.226 

2.95 

Standard 

deviation 

0.99 

0.039 

0.15 

* Under  identical  conditions  and  based  upon  kj,  it  is  approximately  1 5 times  more  reactive  than  paraoxon  and 
1 .6  times  more  reactive  than  amiton.  G.  L.  Horton.  J.  R.  Lowe,  and  C.  N.  Lieskc,  unpublished  results. 
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Spontaneous  reactivation  of  NMPP-inhibited  AChE  involves  displacement  of  the 
phosphinyl  group  with  concomitant  return  of  enzymatic  activity  (equation  2b). 


c 


O 


- P - O - AChE 


CH-i 


Inhibited  AChE 
(inactive) 


spontaneous 
reactivation  , 


Active  AChE  + 


0 

11  r- 
-p-o»" 

1 

CH, 


(2a) 


(2) 

(2b) 


Its  counterpart  has  been  termed  aging  (equation  2a).  With  phosphates  and  phosphonates  this 
pathway  has  been  shown  to  involve  loss  of  an  alkoxy  group  from  the  phosphorus  moiety.7,1 7,18 


In  the  case  of  phosphinylated  AChE  one  should  be  able  to  dismiss  aging  from  reaction 
considerations  (see  INTRODUCTION).  Our  work  to  date  with  NMPP  supports  this  rationale.  A 
summary  of  our  spontaneous  reactivation  results  with  AChE  inhibited  with  NMPP  is  shown  in 
table  3.  They  clearly  reflect  100%  spontaneous  reactivation  at  both  pH  7.60  and  pH  9.10,  with 
t]/2  = 32.2  hours  at  the  lower  pH  and  tj /2  = 9.6  hours  at  the  higher  pH.  By  way  of  contrast, 
the  tj/2  for  the  spontaneous  reactivation  of  phenyl  methylphosphonylated  AChE  under  the 
same  conditions  at  pH  7.60  is  42.5  hours,  with  83%  return  of  activity.5  The  less  than  100%  return 
of  activity  reflects  operation  of  the  parallel  first-order  aging  reaction.  The  trend  of  the  spontaneous 
reactivation  rates  for  the  above-mentioned  compounds  is  in  the  same  direction  as  those  frequently 
found  for  hydrolysis  rates  of  the  analogous  esters;  i.e.,  phosphinates  > phosphonates  > phosphates. 
However,  the  magnitude  of  the  difference  between  these  spontaneous  reactivation  rates  is  not  as  . 
great  as  one  might  expect.  For  example,  p-nitrophenyl  diethylphosphinate  is  14.6  times  more* 
susceptible  to  alkaline  hydrolysis  than  ethyl  p-nitrophenyl  ethylphosphonate.  Likewise,  ethyl 
p-nitrophenyl  ethylphosphonate  hydrolyzes  39  times  faster  in  pH  8.3  buffer  than  the  corresponding 
phosphate,  diethyl  ^-nitrophenyl  phosphate.1  9 These  trends  for  hydrolysis  rates  also  prevail  upon 
replacement  of  the  ethyl  groups  in  the  above-cited  phosphinate  and  phosphate  compounds  by 
n-butyl  groups;  that  is,  at  25.0°C,  the  kQjj  for  the  hydrolysis  of  p-nitrophenyl  c/z-w-butylphosphinate 
is  8 1 .5  M'*  min'* , whereas  the  kQ^j  for  rf/-n-butyl  phosphate  is  0.456  M"*  min'* . This  corresponds 
to  the  phosphinate  being  1 79  times  more  susceptible  to  alkaline  hydrolysis  than  the  analogous 
phosphate. 

Table  3.  Spontaneous  Reactivation  of  Methylphenylphosphinylated 
Eel  Acetylcholinesterase,  25.0°C 


pH 

Replicates 

^obs 

1 1 /2 

Percent 
reactivation  a 

hr'1 

hr 

7.60  b 

4 

0.0215  ± 0.0003 

32.2 

98.9  ± 2.3 

9.10  c 

2 

0.072  ±0.004 

9.6 



103.4  ± 4.0 

a Calculated  using  equation  I. 
b MOPS  buffer. 
c Bicine  buffer. 
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The  statements  just  noted  strongly  suggest  that  the  spontaneous  reactivation  of 
organophosphorus-inhibited  AChE  is  far  more  complex  than  either  general  base  catalysis  or 
nucleophilic  displacement  of  the  phosphorus  moiety  would  predict.  Spontaneous  reactivation 
is  probably  related  to  charge  distribution  around  the  substituents  on  the  phosphorus  atom  in  the 
active  site.  The  electronic  effects  of  substituents  are  influenced  by  the  dielectric  constant  of  the 
reaction  medium.20  Therefore,  one  may  cautiously  speculate  that  a critical  factor  in  spontaneous 
reactivation  is  the  dielectric  constant  in  the  active  site  region  as  influenced  by  the  steric  factors  of 
the  groups  attached  to  phosphorus  in  the  inhibited  enzyme  and  their  perturbation  on  the 
conformation  of  the  enzyme. 

Our  oxime  reactivation  experiments  confirmed  the  stoichiometry  obtained  in  our 
spontaneous  reactivation  studies  and  they  support  the  absence  of  aging  in  NMPP-inhibited  enzyme. 
Both  TMB-4  and  2-PAM  completely  reactivated  methylphenylphosphinylated  enzyme  in  1 hour. 

As  expected,  the  same  concentration  of  MINA  was  found  to  be  less  efficacious;2 1 that  is,  MINA 
has  been  reported  to  be  only  about  1/7 Sth  as  effective  as  2-PAM  in  the  reactivation  of  diethyl- 
phosphoryl  AChE  (human  erythrocyte).  The  structural  formulas  of  these  oximes  and  the  results 
of  these  investigations  are  shown  in  table  4. 

Table  4.  Oxime-Induced  Reactivation  of  NMPP-inhibited  Eel  AChE, 

[Oxime]  = 1 X IQ"4  M,  pH  7.60,  0.10  M MOPS  Buffer,  25.0°C 


OXIME 

STRUCTURAL  FORMULA 

REPLICATES 

% REACTIVATION 
IN  1.0  HOUR 

2-PAM 

P) 

ch3 

L" 

L_  c = N — OH 
Cl© 

2 

96.4  + 0.2 

TMB-4 

Cl© 

H 

1 

c = 

0 

7© 

ch2 

H 

1 

M — OH  C = N — OH 
1 

jP*| 

7®  Cl© 

— ch2  — ch2 

2 

99.8  + 0.8 

MINA 

O H 

ii  1 

CH3  — C — C=N  — OH 

2 

28.5  + 0.4 
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Since  induced  reactivation  by  MINA  was  incomplete  at  the  end  of  1 hour,  the  reaction 
was  monitored  through  a 4-hour  period  in  two  experiments.  The  apparent  second-order  reactivation 
rate  constant,  k2(app)’  was  calculated  using  equation  3 

^obsd 

k2(app)  = [OXIME] 

where  kobsd  is  the  observed  first-order  reactivation  in  the  presence  of  an  oxime  concentration.  The 
value  obtained  was  48.2  M'1  min'*  ± 0.51%,  no  distinction  being  made  between  protonated  and 
unprotonated  forms  of  the  oxime.  Likewise,  cognizance  is  taken  that  the  reaction  is  not  truly 
second  order  but  that  it  involves  an  intermediate  complex  of  the  oxime  with  the  inhibited  enzyme. 

Thus,  although  k2(appj  values  give  only  a rough  measure  of  the  reactivating  effectiveness 
of  a compound,  it  frequently  suffices  for  qualitative  screening  purposes.22  As  such,  it  is  interesting 
to  note  that  the  MINA-induced  reactivation  of  eel  AChE  inhibited  with  NMPP  is  more  effective 
than  MINA  reactivation  of  human  erythrocyte  cholinesterase  inhibited  with  Sarin  (21 .5  M*1  min'1 ) 
or  tetraethyl  pyrophosphate  (7.2  M'1  min*1).23 

In  conclusion,  this  report  documents  (l)methylphenylphosphinylatedeel  AChE  as  a model 
system  suitable  for  use  in  examining  the  effects  of  various  perturbing  agents  on  the  spontaneous 
reactivation  process  and  (2)  the  first  spontaneous  reactivation  data  on  a phosphinate-inhibited 
cholinesterase.  We  are  currently  in  the  process  of  examining  the  effects  of  various  compounds  on 
the  spontaneous  reactivation  of  methylphenylphosphinylated  eel  AChE.  In  addition,  several  additional 
phosphinate  esters  are  being  synthesized,  which,  based  upon  linear  free  energy  relationships,2  4 will 
display  reactivating  characteristics  different  from  those  found  in  this  study.  The  use  of  such 
compounds  should  allow  greater  flexibility  in  spontaneous  reactivation  studies  of  inhibited  AChE 
and  ultimately  contribute  to  the  therapy  of  organophosphate  poisoning. 
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APPENDIX  A 


DOC  UMENTATION  FOR  PROGRAM  EXPFIT 

A.l.  Macrodocumentation  of  Program  EXPFIT 

a.  Title:  MEDICA*BIOCHEM. EXPFIT 

b.  Programmer : CPT  John  R.  Lowe,  Biomedical  Laboratory 

(SAREA-BL-RE) , 671-3836  or  2626. 

c.  Machine  and  Language:  UNIVAC  1108,  FORTRAN  V. 

d.  Purpose:  This  program,  a modification  of  one  written  by  W.L.  Wilcoxsonl 

will  calculate  by  an  iterative  least  squares  method  the  observed  first-order 

rate  constant  and  the  limiting  velocity  for  the  spontaneous  reactivation  of 
inhibited  enzyme  using  equation  1 and  also  will  generate  a computer  graphics 
datafile  for  subsequent  off-line  automated  plotting.  This  program  will  accom- 
modate up  to  200  sets  of  velocity  versus  time  data  for  each  rate  constant  to 
be  calculated.  A flowchart  for  this  program  is  given  in  section  A. 2.  A copy 
of  the  program  is  listed  in  section  A. 3. 


Input  Parameters: 

IDRUN 

experiment  identification 

ID2 

continuation  of  IDRUN 

ID3 

continuation  of  IDRUN  and  ID2 

I BLANK 

dummy  variable  to  accept  a value  in  the  datafile 

which  is  used  in  an  alternate 

computer  program 

IB2 

continuation  of  IBLANK 

IB3 

continuation  of  IBLANK  and  IB2 

TUN  IT 

units  of  time  (hours,  minutes. 

or  seconds) 

VUNIT 

velocity  function  (%  reactivation,  absorbance, 
% activity,  or  velocity) 

VU2 

continuation  of  VUNIT 

VU3 

continuation  of  VUNIT  and  VU2 

TIME(I) 

time  elapsed  at  which  velocity 

function  was  measured 

RATE (I) 

velocity  function  measured  at 

time^ 

Input  Format 

Line/Card  No 

. Parameter  Columns 

Format  Specification 

1 

IDRUN  1-4 

A4 

1 

ID2  5-8 

A4 

1 

ID3  9-12 

A4 
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Line/Card  No. 

Parameter 

Columns 

Format  Specification 

2 

TUN  IT 

1-6 

A6 

3 

VUNIT 

1-6 

A6 

3 

VU2 

7-12 

A6 

3 

VU3 

13-18 

A6 

4 

TIME(l) 

1-8 

F8.3 

4 

RATE(l) 

9-17 

F9.4 

5 

TIME(2) 

1-8 

F8.3 

5 

RATE(2) 

9-17 

F9.4 

N+4 

TIME (N+l) =0.0 

1-8 

F8.3 

N+4 

RATE (N+l) =0.0 

9-17 

F9.4 

N+5 

IDRUN='N0NE' 

1-4 

A4 

g.  Output  Parameters:  In  addition  to  each  of  the  parameters  listed 

in  section  A.l.e,  the  following  parameters  are  printed  out  in  one  of  two 
output  files,  a hard-copy  analytical  summary  (unit  6)  and  a plotting  data- 
file (unit  16). 

DELTA  predetermined  small  number  used  to  terminate  iterations 

based  on  evaluating  percentage  change  between  successive 

estimations  of  k , , V„,  and  V 

obs  0 °° 

HALF  calculated  half  life  of  the  reaction 

IT  iteration  number 

ITMAX  maximum  number  of  iterations  allowed 

M calculated  plotter  control  digit  for  plotting  program 

N number  of  sets  of  time  versus  velocity  data 

NN  number  of  plotting  points  generated  for  plotting  datafile 

OBSK  calculated  observed  rate  constant 

QHALF  number  of  half  lives  over  which  the  data  were  collected 

SAS  analytic  variance  for  the  fitted  curve 

SS  experimental  variance  for  the  fitted  curve 

SOD  sum  of  differences  between  experimental  and  calculated 

velocities 

TCALC  time  calculated  by  fractionation  of  total  experimental 

time  interval,  used  for  generating  plotting  data 

TIMO  abscissa  intercept,  time  when  velocity  = 0 
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VDIFF  difference  between  experimental  velocity  and  the 

velocity  calculated  from  the  fitted  equation 

VCOMP  velocity  calculated  from  the  fitted  equation 

VINF  calculated  maximum  velocity  at  infinite  time 

VZERO  ordinate  intercept,  velocity  at  time  = 0 

h.  Output  Format: 

(1)  Unit  6 — Hard  Copy  Analytical  Summary 

A copy  of  this  output  is  shown  in  section  A. 4 

(2)  Unit  16 — Plotting  Datafile 


Line/Card  No. 

Parameter 

Columns 

Format  Specification 

1 

IDRUN 

1-4 

A4 

1 

ID2 

5-8 

A4 

1 

ID  3 

9-12 

A4 

2 

TUNIT 

1-6 

A6 

3 

VUNIT 

1-6 

A6 

4 

OBSK 

1-10 

E10.4 

5 

TIMO 

1-10 

E10.4 

6 

TIMO 

1-15 

E15.8 

7 

VINF 

1-10 

E10.4 

8 

VINF 

1-10 

E10.4 

9 

HALF 

1-15 

E15.8 

10 

NN 

1-3 

13 

11 

TIMO 

1-8 

F8.3 

11 

VCOMP 

9-17 

F9.4 

12 

TCALC 

1-8 

F8.3 

12 

VCOMP 

9-17 

F9.4 

210 

TCALC 

1-8 

F8.3 

210 

VCOMP 

9-17 

F9.4 

211 

M 

1-3 

13 

212 

N 

1-3 

13 

213 

TIME(l) 

1-8 

F8.3 

213 

RATE(l) 

9-17 

F9.4 

• 

N+212 

TIME(N) 

1-8 

• 

• 

F8.3 

N+212 

RATE(N) 

9-17 

F9.4 

N+213 

M 

1-3 

13 

N+214 

NN 

1-3 

13 

N+215 

TIMO 

1-8 

F8.3 

N+215 

VCOMP 

9-17 

F9.4 
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Line/Card  No. 

Parameter 

Columns 

Format  Specification 

N+216 

TCALC 

1-8 

F8.3 

N+216 

VCOMP 

9-17 

F9.4 

N+414 

TCALC 

1-8 

F8.3 

N+414 

VCOMP 

9-17 

F9.4 

N+415 

M 

1-3 

13 

N+416 

N 

1-3 

13 

N+417 

TIME(l) 

1-8 

F8.3 

N+417 

RATE(l) 

9-17 

F9.4 

2N+416 

TIME(N) 

1-8 

F8.3 

2N+416 

RATE(N) 

9-17 

F9.4 

2N+417 

M 

1-3 

13 

2N+418 

IDRUN= 'NONE ' 

1-4 

A4 

h.  Operating  Instructions: 

(1)  For  batch  operations,  follow  the  instructions  listed  in  sections 
3,  7,  16,  and  17  of  the  Edgewood  Arsenal  UNIVAC  1108  User's  Guide. 

(2)  For  demand  operations,  the  required  UNIVAC  1108  system  inter- 
actions are  described  in  the  following  section.  All  of  our  data  reduction  is 
performed  in  the  demand  mode.  The  telecommunication  port  to  the  computer  is 
initialized  by  calling  the  appropriate  phone  number  for  the  baud  rate  being 
used  and  by  typing  in  the  site  identification  number.  After  initialization, 
the  demand  runstream  is  begun  with  an  Executive  8 control  statement  (master 
space  RUN).  This  control  statement  contains  the  run  identifier,  account  number, 
project  identifier,  and  time/page  estimates.  Experimental  data  are  then  entered 
into  an  element  of  an  on-line  datafile  (MEDICA*EXPDATA. ) in  the  format  listed 

in  section  A.l.f  by  using  the  ELT  processor  for  program  file  maintenance. 
Proofreading  and  subsequent  corrections  to  this  data  element  are  made  by  using 
the  editing  processor  (@ED,U  FILE. ELEMENT) . After  proper  format  has  been 
established  for  the  data  element,  the  data  reduction  is  begun  by  typing  in  the 
following  set  of  system  control  commands. 

@ERS 

@ASG,T  TEMP. ,F2 

@USE  16.  TEMP. 

@FOR,N  MEDICA*BIOCHEM. EXPFIT ,TPF$ . 

@MAP,N 

IN  TPF$ . 

@XQT 

In  actual  practice,  these  statements  have  been  stored  as  a canned  runstream 
in  the  file  MEDICA*RUNEXPFIT . Therefore,  in  order  to  start  the  data 
processing  by  the  FORTRAN  curvefitting  program,  we  type  in  @ADD  MEDICA*RUNEXPFIT. 
Now  the  only  requirement  to  complete  the  data  reduction  is  to  answer  the 
questions  which  are  printed  out  at  the  teletype.  A sample  of  the  output  is 
shown  in  section  A. 4. 
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A 2.  Flowchart  tor  Program  EXPFIT 
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A. 3.  FORTRAN  Statements  for  Program  EXPFIT 


C....  DIC4*8ICCHEH.  EXPFIT  

C MODEL  V = V INF  • Cl  - EXPI-K  • IT  - TO))3 


130  FORMAT  1/OXtTl.'  .*•••••••••••••••••••••••••••••••••••»•••••••••• 

••••••••,/  ) 

110  FORMAT  I 3 X. T1 « * IF  YOUR  DATA  HAS  BEEN  PREVIOUSLY  STORED  IN  AN  ON-L 
TINE  FILE,*/*  ENTER  »9AD0  F ILEN  AME.  ELEMEN  INANE*  . FOR  EXAMPLE:  8A  Of) 
aEXPOATA.f X2854  */ • OTHERWISE.  TYPE  IN  YOUR  RUN  IDENTIFICATION.*) 

123  cCRMAT  (JAM 

133  FORMAT  ( 3 X ♦ T 1 « • WHAT  ARE  THE  UNITS  OF  TIME  I HRS  # MIN,  OR  SEC)?*) 
143  FORMAT  < A 6 ) 

153  FORMAT  (A3) 

163  FORMAT  1 3 X » T1 » • WHAT  IS  THE  VELOCITY  FUNCTION  I*  REACTIVATION,  A3S 
S0P8ANCE,  X ACTIVITY, OR  VELOCITY)?*) 

173  FORMAT  (3A6) 

183  FORMAT  (3X.T1.*  ENTER  TIME  * VELOCITY  IN  THE  FOLLOWING  FORMAT:*/* 
ATTTT. TTT VVVV.VVVV*  ) 

193  FORMAT  (F9.3.F9.4I 

203  FORMAT  (1H3) 

213  FORMAT  (8HDDELTA  : »E  1 0 . 4 . 5 X . 7H  IT  MA  X : . 14 , / . 1 3H OC YCLE  V(INF)./ 

%*  ♦*♦♦♦ 

223  FORMAT  1 1 X , I 3 , E 15. 5 ) 

2 33  FORMAT  ( 6 5HQT  HE  PROGRAM  HAS  CYCLED  THROUGH  THE  MAXIMUM  NUM3ER  OF  I 
ITERATIONS.) 

243  FORMAT  C10H3RUN-ID  : »2X.3A4,/»9H  KI09S)  : . E 1 S .8 » 1 OX  ,6HT  ID)  :,E15. 
4 6 • / • 9H  VI  INF)  :*E15.8«1QX.6HV(  0 A :.E15,8./.9H  TII/2)  = . E 1 S . 8 » 1 X . A 6 
W./9H39ASE0  0N.I3.16H  POINTS  CO  VE  RI  NG  .E  15 .6 . 1 2H  HALF  LIVES.) 

253  FORMAT  IE13.4) 

263  FORMAT  IE15.8) 

2 73  FORMAT  (1H  .1X.7HTIME  IN . 1 3X  . 1 2hE X P£ RI ME  NT AL . 8 X. 10HC AL CULA TED . 8 X . 

413hDIFFFRENCE,/T1»4X»A6#9X,3A6»2X.3A6,  10 HI E X P-CALC ) • /Tl, * ♦♦♦♦♦♦♦♦ 
«♦•. 13X. T2 3. •♦♦♦♦♦♦♦♦♦♦♦♦♦♦*, 7X.T40.* ♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦».SX.T 58 •*♦♦♦♦♦♦ 
?♦♦♦♦•  ) 

283  FORMAT  IT1»F8.3*T2Q»F8.3,T39»F8.3,T57,F8.3) 

293  FORMAT  ( 24H0EXP£RI MENTAL  VARIANCE  :»E15.8./24H  ANALYTIC  V AR I AN 

* CE  :.E15.8,/24H  SUM  OF  DIFFERENCES  z.ElS.8) 

330  FORMAT  ( 13) 

313  FORMAT  I3X.T1,*  NEXT  RUN-ID?*) 

323  FORMAT  I3X.T1,*  AN  EOUATION  HAS  BEEN  FITTED  TO  YOUR  DATA.*/*  NOV  T 
4Y»E  IN  -8C0PY.I  TE MP. • MEO IC A *E  XPF I T*LOT. *, 3A4, *■ •) 

DIMENSION  TIMEI200),RATEC230) 

DELTA:. 0303 1 
I TM A X :4  0 
WRITE  16.130) 

WRITE  16,110) 

3 30  READ  15.123)  I ORUN . I 02 , ID  3 

WRITE  116.120)  IDRUN.I02.IC3 
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O O O Cl  Cl 


340  IF  ( I0RUN.E0. *N0NE *t  GO  TO  540 
IP4=IPRUN 
105=102 
106=  103 
N =0 

WRITE  16*130) 

REAO  (5*140)  TUN  IT 
WRITE  (16*150)  TUNIT 
WRITE  (6*160) 

READ  (5*170)  VUN IT • VU2*  VU  3 
WRITE  (16*170)  VUN IT 
WRITE  (6*180) 

00  350  1=1*200 

READ  (5*193)  T IHE ( I) . R AT E ( I ) 

IF  (RATE( I) .EQ.O.)  GO  TO  360 
N ;N*  1 

350  CONTINUE 
360  WRITE  (6*200) 

WRITE  (6*210)  OELTA*ITMAX 
V2ER0=RATE(  1) 

TEST  = .5*(RATE(N)*VZER0  ) 

CO  370  1=3* N 

IF  <RATE( I) .GT.TEST)  GO  TO  380 
370  CONTINUE 

03SK0=1./TIHE(N) 

GO  TO  390 

380  03SK  C = 1 . / TIME ( I ) 

390  VINFO=RATE(N) 

IT  = 0 

40C  CU=0.0 
C 1 2 =0 . 0 
C 1 3 =0 . 0 
C 1 R =0 .3 
C22  =3 . 3 
C23=0.0 
C2R  =3 . 0 
C3  3 =0 . 0 
C3R  =3 . 0 
CRR=0.3 
DC  430  1=1. N 


EXPON=-OBSKC*TIME( I) 

IF  (EXP0N.GT.-49.)  GO  TO  410 
U = 1.0 
V =0 .0 
W =0 .3 
GO  TO  423 
413  V=EXP(EXPON) 

U=1 .- V 
W=TIPE(I)*V 
423  TEST=RATE(I> 


I 
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cn=cn*u*u 

C12=C12*U*V 
C13  = : 13*U»W 
c 1 r=c ir*u*test 
C22=C22*V*V 
C23=C23»V*W 
C2R=C2R*V*TEST 
C3  3 33  ♦ W • V 

C3Rrr  3R*W*TEST 
CRR=CRR*TEST*TEST 
4 30  CONTINUE 

A;C11*C33-C13»C13 
3=C12*C33-C13*C23 
C-C22*C33-C23»C23 
E=C1R»C33-C3R*C13 
E=C2R*C  33-C3R*C23 
VINF=(E*C-E»B)/(A*C-B*e) 

V7ER=(  F-VINF*  B) /C 

Q=(C3R-VINF*C13-VZER*CZ3)/C33 

0ESK=03SKO*O/(VINfO-V?ERO) 

I T=IT*1 

OKTEST:|. -OBSK/OBSKO 
03SK0=08SK 

V TEST  =1 .-V2ER/VZER0 

V 7E  RC - V ZF  p 

VITEST=1.-VINF/VINF0 

VINE0=VINE 

IF  (DELTA. LE.  ABS(OKTE«>T)  ) GO  TO  440 
IF  (DELTA. LE. ABS(VTEST  )>  GO  TO  440 
IF  (DELTA. LE. ABS(VITEST)  ) GO  TO  440 
50  TO  453 

443  WRITE  (6*  220)  IT.VINF 

IF  ( IT.LT.ITMAXJ  GO  TO  433 
WRITE  <6.2331 

4 53  TIM3=( ALOG(  l.-VZER/VINE)  )/OBSK 
HALF=.S9315/09SK 
3HALF=TIME(N) /HALF 

WRITE  ( 6. 2401  IDRUN.  102 . 103  . OBSK  .T  INO.  VI  NF  tVZER.HALF  .TUN  IT  .N.QKALF 
WRITE  (6.230) 

WRITE  <16. 253)  03SK 
WRITE  (16.253)  TIMO 
WRITE  (16.263)  TIMO 
WRITE  (16.253)  VINE 
WRITE  (16.253)  VINE 
WRITE  (16.263)  HALF 

WRITE  (6.270)  TUNIT.VUNIT.VU2.VU3.VUNIT. VU2.VU3 
ENM3rFLOAT(M-3> 

SAS=( CRR-VINF*CIR- VZEP*C2R-0*C3R )/ENM3 

50S=3.3 

SC0=3 .3 
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<♦60 
4 70 

4 5 0 

4 S3 

533 

513 

523 

5 30 
540 


■■ 


— 


00  470  1 = 1. N 
F XP0N=-09SK*T IME  <1  ) 

VCCMPrVINF 

IF  (EXP0N.LE.-49.)  GO  TO  460 
VCOMP=VCOMP-<  VINF-VZEB )*EXP1FX  PON) 
VOIFF=PATE( I) -VCOMP 
SCS=SOS*VDIFF*VOIFF 
SOD=SOO*VOIFF 

WRITE  (6»2e0>  TIME  (I  ). RATE  (I  »•  VCOMP,  VO  IFF 

CONTINUE 

5S=SCS/ENM3 

WRITE  (6.290)  SS.SAS.SOO 
T I MN  = T IMF ( N ) 

T I NF  = 5.  • HALF 
NN=233 

WRITE  (16.303)  NN 
VCOM’rO .3 

WRITE  (16.193)  TIMO, VCOMP 
DO  493  1=1.199 

TCALC=TIMO*FLOAT (I )* (T IMN-TIMO  )/ 199. 

VCGMP  = VINE* (1 .-EXP (-09SK*( TCALC-TIMO  ))  ) 

WRITE  (16.193)  TCALC. VCOMP 

CONTINUE 

m=0 

WRITE  (16.303)  M 
WRITE  (16.303)  N 
DC  50  3 J=  1 » N 

V ° I T E (16.193)  TIME(J).RAT£(J) 

CONTINUE 
M = -N/  33 

IF  <*>  520.  51  0.  51  3 
M = -l 

WRITE  (16.303)  M 

IF  ( TIMN.rO. TINF  ) GO  TO  5 30 

T I MN  = T I Nf 

GO  TO  483 

WRITE  (6.130) 

WRITE  (6.310) 

RE  AC  (5.1231  iaLANK.I3L2.I3L3 
IF  ( I9LANK.E9. * •)  GO  TO  330 

I!3RUN  = I3L  AN* 

IC2=I3L2 
1 03  = I 3L  3 
GO  TO  343 
WRITE  (6.130) 

WRITE  (6.323)  ID4.I05.I0E 

WRITE  (6.133) 

3TCR 

ENC 
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A.  4.  Sample  of  Teletype  Interactive  Output  for  Program  EXPF1T 

a.  The  following  Is  a sample  of  what  Is  seen  during  a demand  operation  at  the 
teletype  terminal.  The  first  input  Is  the  site  Identification.  The  system  responds 
with  the  Edgewood  Arsenal  identification  line.  The  operator  then  enters  his  run  control 
command.  The  system  responds  with  the  date  and  time  and  then  indicates  that  it  is  ready 
for  input  with  a ">"  sign  in  column  one.  Each  time  a ">"  sign  is  printed,  the  operator 
enters  his  control  commands  or  data. 


: ? I T E~ I j 

:*EDGEMOQD  ARSENAL  1108  SYSTEM  "FF:.  31-844  IJPD  HIRSH* 
0F  i IH  .JRLOME > COST-CODE  - PROJECT- 1 D > 30 » 1 000 
DftTE:  088178  TIME:  1 85488 

>©ADIi  MED  I CA*RUNEXF‘F  I T . 

FlJPPUP  M088-08--  8 1-18: 54 

READY 

READY 

Fi  .p  : .0E3-0t.-  '81  • '78-18: 54:  44  (0-  ) 

END  FOP 

MAP28P1  F‘L71-3  08-"81.'?6  16:54:58 
END  NAP 


IF  • ni.lp  DflTfi  HRS  DEEM  PPEUIOUSLY  STORED  IN  RN  ON-LINE  FILE- 
ENTER  "MR DIi  F I LE NOME . ELEMENTNRME " . FOR  EXAMPLE:  ©ADD  EHF'DATR.  EX8854 
OTHERWISE-  TYPE  IN  YOUR  RUN  IDENTIFICATION. 

>C?ADD  EJ-IPDATR.  EH8854 

NHRT  APE  THE  UNITS  OF  TIME  (HRS-  MIN,  OR  SEC)? 

HURT  IS  THE  UELOCITY  FUNCTION  (X  REACT I MAT I ON - ABSORBANCE-  % ACTIMITY- 
OP  "ELOCITYr? 

IT  E S' F TIME  f ML'LOCITY  IN  THE  FOLLOWING  FORMAT: 

TT7  T , TTTU1  n h i.  1 11 11  in 


DF.L  A 


i 000-04 


1'  TIT 


1 1 ( It  ID. 

33804+08 
99985+08 
1 0007+03 
i -'11107+03 
1 t'i007+03 


itna;-:  ~ 


rum- id  - 

MOBS 'I  - 
It!  IMF  I - 

Tll/2)  - 

: tiSE’.  ufi 


L:  :P8854 
. 3 1 1 99663-0 1 
, 1 000^400+03 
346968 74+08  HITS 

9 POINTS  COHERING 


40 


T f CO  » -.85671^68+00 

11  ( 0 ;i  :=  .1 80 1 1 5O9+0 1 

.308139+01  HALF  LIMES. 
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r "—i 


’ I ME  IM 

EXPERIMENTAL 

CALCULATED 

DIFFERENCE 

HRS 

•.  REACT I NATION 

’.  REACT I NATION 

i EXP-CALC ) 

♦♦  ♦+♦♦+ 

++♦♦♦♦♦++♦++♦+ 

■*■++++♦■*•+♦+++++ 

+++♦+ + 

1 . O00 

3.05O 

3.863 

• '•  * 

".04  2* 

1 . 080 

3. 560 

4 . 026 

- . 466 

2 . 1 60 

5.940 

8.200 

- . 260 

> . 000 

7.640 

7.357 

-.217 

4 . 050 

10. ISO 

9 . 88  7 

. 293 

5.  160 

12.050 

1 1 ! 984 

. U66 

6 . 000 

13.240 

13.539 

- . 299 

" . 000 

15. 280 

15.354 

-.074 

,.250 

1 1 • 820 

17.570 

. 250 

i 1 . 050 

22. 1O0 

L.L.  • i 1 

-.225 

12.1 60 

24.650 

24. 133 

.517 

13.080 

26. 7O0 

25 . 600 

1 . 100 

1 4 . 050 

27. 330 

27.  115 

. 265 

1 5 . 5O0 

30. 1O0 

29. 324 

. 77 6 

2 1 . 500 

36. 340 

87. 77 4 

- 1 . 484 

23. 000 

41.330 

39.724 

1 . 606 

24 . 580 

42 . 680 

41.712 

. 968 

26 . 500 

44 . 1 30 

44.040 

. 090 

2-. 660 

44.550 

45.401 

-.351 

23 . 5O0 

46.540 

46 . 366 

. 174 

29.500 

46. 710 

47. 493 

— "*83 

3 1 . 000 

43. 350 

49. 1 39 

- . 7 89 

32 . -'6*  1 

50. 340 

50 . 586 

-.246 

33. 060 

50 . 560 

51.315 

-.755 

36.330 

54.450 

54.531 

-.  131 

38. 910 

59.470 

57. 002 

2. 468 

45.50O 

61.1 70 

62 . 6 1 8 

-1.448 

4?. 420 

65. 120 

64 .112 

1 . O08 

52.250 

69. 480 

67. 612 

1 . 868 

55 . 660 

73. 820 

69.876 

3.  944 

70. 160 

73. O30 

77 .368 

-4.338 

70 . 6 1 0 

76.540 

TO  £4  TO 

-1.538 

72. 330 

75. 09O 

78, 866 

-3. 776 

72 . 660 

76. ISO 

79.014 

-2.334 

76 . 950 

82. 110 

80.845 

1 . 265 

30 . 000 

81 . 480 

32. 049 

- . 569 

94 . 920 

88. 960 

86.936 

2.024 

97 . 0O0 

90. 370 

87 . 503 

c • ob  i 

100. 750 

89. 260 

88.464 

. 796 

EXPERIMENTAL 

UARIAMCE  = .23565553+01 

ANALYTIC 

UARIANCE  = .28565209+01 

SUM  OF  DIFFERENCES  = .33676624-04 

».3«i«9«IK3«3«)K3KIK9KI«JK9KltHK3«3«3«MC3«3tHKI«n«9ti»gi«3«3«3«)K3K3«l}KI«iK»H«3«IJKiKiKilHK3K«3liM(i«linik»lKMMWI 


next  rum- id? 

NOME 


AN  equation  has  deem  fitted  to  your  data. 

MOW  TYPE  IM  "0COPY*  I TEMP . » MED  I CA*EXPF  I TPLOT . EXP2S54 

******KWKK*K*ft********K***K*«********M«*K*WjK*KSK*KKKKSK!i;WSMM**MMKKK!KIKKM 


>0COPY , I TEMP . > MED I CA»EXPF I TPLOT . EXP2854 
FUPPUR  O026-06/2 1-16:56 
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b.  In  the  previous  example,  the  data  were  entered  by  using  an  on-line  datafile.  If 
data  are  to  be  entered  during  program  operation  manually,  the  following  is  a sample  of  the 
portion  of  the  runstream  which  is  different  from  that  shown  in  A.U.a.  Instead  of  responding 
to  the  first  statement  with  §ADD  EXPDATA.EX285** , the  operator  would  type  in  the  run  identi- 
fication, EXP285U.  The  other  questions  are  then  answered  in  turn  as  the  ">"  sign  appears. 


E:  :P£854 

UHAT  ARE  THE  UNITS  OF  TIME  (HRS*  MIN*  OR  SEC)? 

>HRS 

WHAT  IS  THE  "ELGCITY  FUNCTION  C.  REACT I NATION*  ABSORBANCE*  ACTIUITY* 
OR  "ELOCITYJ? 

>*.  REACT  I MAT  I ON 

ENTER  TIME  & "ELOCITY  IN  THE  FOLLOWING  FORMAT: 

TTTT.TTT"'  'UU.UU'H' 


1 . 00 

3.05 

/' 

1 . 08 

3 . 56 

\ 

c!.  lb 

5.94 

3.00 

7.64 

■> 

4 . 05 

10.  18 

V. 

5.  IS 

1 £ . 05 

^•5. 

6.  OO 

13.  £4 

7.00 

15.  c!8 

8 . £5 

1 7 . 3£ 

i 1 . 05 

1 3c.  1 0 

y 

1 £ • It 

£4.65 

13.  Oy 

£6 . 70 

14.05 

£7.33 

> 

15.50 

30.  10 

> 

£ 1 . 50 

36 . 34 

> 

£3.00 

41.33 

£4.58 

4£ . 68 

y 

£t..50 

44.  13 

> 

£7.66 

44.55 

]> 

£3.50 

46.54 

> 

£9.50 

46.  71 

"S 

31.00 

48.  35 

> 

3£  ■ 36 

50  ■ -34 

y 

33. 06 

50.56 

y 

36. 33 

54.45 

38.  91 

59.  47 

V 

45.50 

61.17 

47.  4£ 

65.  1£ 

\ 

5£.  £5 

69.43 

/ 

55. 66 

73. 8£ 

> 

70.  16 

73.03 

y 

70.61 

76. 54 

y 

1 t-  • 

75.09 

y 

72.  66 

76. 18 

y 

76.95 

S£.  11 

"y 

80.00 

81.48 

94.  9£ 

83.96 

/ 

97.00 

90.  37 

y 

100.75 

89.  £6 

' 

0.0 

0.  0 
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DOCUMENTATION  FOR  PROGRAM  EXPF1TPLOT 

B.l.  Macrodocumentation  for  Program  EXPFITPLOT 

a.  Title:  MED I C A* B1 OCR EM. EXPFITPLOT 

b.  Programmer:  CPT  John  R.  Lowe,  Biomedical  Laboratory 

(SAREA-BL-RE) , 671-2626/3836. 

c.  Machine  and  Language:  UNIVAC  1108,  FORTRAN  V. 

d.  Purpose : This  program  will  provide  the  required  system  control 

commands  to  plot  the  calculated  least-squares-fitted  exponential  curve 
from  program  EXPFIT.  The  off-line  CALCOMP  plotter  at  the  main  computer 
site  and  the  library  plotting  subroutines  for  this  plotter  are  used  in  this 
program.  Up  to  200  sets  of  abscissa-ordinate  data  can  be  accommodated 

to  draw  smooth  curves  for  the  fitted  exponential  equation.  A copy  of  the 
flowchart  for  this  program  is  shown  in  section  B.2.  A copy  of  the  program 
is  listed  in  section  B.3. 

e.  Input  Parameters : 

HALF  calculated  half  life  of  the  reaction 


IDRUN 


experiment  identification 


continuation  of  IDRUN 

plotter  control  digit  used  to  selectively  choose 
the  number  of  points  to  draw  on  the  plot 

number  of  sets  of  TIME  versus  RATE  data 

observed  rate  constant  in  alphanumeric  format 


OBSKA 


continuation  of  OBSK 


OBSKB 


continuation  of  OBSK  and  OBSKA 


RATE(I) 


velocity  function  measured  at  time. 


TIME(I)  time  elapsed  at  which  velocity  function  was 

measured 

TIMO  abscissa  intercept,  time  when  velocity  * 0, 

in  alphanumeric  format 


TIMOA 


continuation  of  TIMO 


TIMOB 


continuation  of  TIMO  and  TIMOA 


TIMO  in  numeric  format 
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TUNIT 

units  of  time  (hours,  minutes, 

or  seconds) 

VINF 

calculated  maximum  velocity  at 

infinite  time  in 

alphanumeric 

format 

VINFA 

continuation 

of  VINF 

VINFB 

continuation 

of  VINF  and  VINFA 

VINFN 

VINF  in  numeric  format 

f.  Input  Format: 

Line/Card  No. 

Parameter 

Columns 

Format 

1 

IDRUN 

1-6 

A6 

1 

ID2 

7-12 

A6 

2 

TUNIT 

1-3 

A3 

3 

VUNIT 

1-6 

A6 

4 

OBSK 

1-3 

A3 

4 

OBSKA 

4-7 

A4 

4 

OBSKB 

8-11 

A4 

5 

TIMO 

1-3 

A3 

5 

TIMOA 

4-7 

A4 

5 

TIMOB 

8-11 

A4 

6 

TIMON 

1-15 

E15.8 

7 

VINF 

1-3 

A3 

7 

VINFA 

4-7 

A4 

7 

VINFB 

8-11 

A4 

8 

VINFN 

1-10 

E10.4 

9 

HALF 

1-15 

E15.8 

10 

N-200 

1-3 

13 

11 

TIME(l) 

1-8 

F8.3 

11 

RATE(l) 

9-17 

F9.4 

12 

TIME(2) 

1-8 

F8.3 

12 

RATE (2) 

9-17 

F9.4 

210 

TIME(N) 

1-8 

F8.3 

210 

RATE(N) 

9-17 

F9.4 

211 

M 

1-3 

13 

212 

N 

1-3 

• 

13 

N+212 

TIME(N) 

• 

1-8 

F8.3 

N+212 

RATE(N) 

9-17 

F9.4 

N+213 

M 

1-3 

13 

N+214 

N=200 

1-3 

13 

N+215 

TIME(l) 

1-8 

F8.3 

N+215 

RATE(l) 

9-17 

F9.4 
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Line/Card  No.  Parameter 


Columns 


Format  Specification 


N+414 

TIME(N) 

1-8 

F8.3 

N+414 

RATE(N) 

9-17 

F9.4 

N+415 

M 

1-3 

13 

N+416 

N 

1-3 

13 

N+417 

TIME(l) 

1-8 

F8.3 

N+417 

RATE(l) 

9-17 

F9.4 

2N+416 

TIME(N) 

1-8 

F8.3 

2N+416 

RATE(N) 

9-17 

F9.4 

2N+417 

M 

1-3 

13 

2N+418 

IDRUN='N0NE 

' 1-6 

A6 

g.  Output  Parameters:  All  output  from  this  program  is  through  the  library 

subroutines  available  on-line  to  provide  CALCOMP  plotter  control.  In  addition 
to  the  parameters  listed  in  B.l.e,  the  following  parameters  are  also  printed. 

ASLEN  length  in  inches  of  the  asymptote  line  at  VINFN 

ASYMP  length  in  inches  to  the  ordinate  value  of  VINFN 

ZERO  length  in  inches  of  time  * 0 on  the  abscissa  axis 

PLOTID  character  used  to  identify  the  two  graphs  described 

h.  Output  Format:  Section  B.4.  contains  the  graphical  CALCOMP  output  which 

results  from  the  operation  of  this  program.  For  each  fitted  exponential  equation 
two  plots  are  drawn.  The  first  has  an  axis  scale  based  on  the  range  of  time 
values  from  TIMON  to  TIME(N)  and  the  second  has  an  axis  scale  based  on  the  range 
in  velocity  values  from  zero  to  VINFN  and  the  range  of  time  values  from  TIMON 

to  five  times  the  calculated  half  life  for  the  reaction.  The  sample  graphs 
show  the  differences  in  graphical  output. 


i.  Operating  Instructions: 

(1)  For  batch  operations,  follow  the  instructions  listed  in  the 
UNIVAC  1108  User's  Guide. 

(2)  For  demand  operations  the  required  UNIVAC  1108  system  interactions 
include  site  and  runstream  initialization  as  described  in  appendix  A,  section 

A. 1.1. (2).  To  initiate  plotting,  a canned  runstream  located  in  MEDICA*PLOfEXPFIT 
is  entered  at  the  teletype  keyboard.  This  file  contains  the  following: 

@ERS 

@ASG,T  TEMP . ,F2 
@USE  16., TEMP. 

@F0R,N  MEDICA*BI0CHEM. PLOTTEREXPFT ,TPF$ . 

@MAP ,N 
IN  TPF$ . 

@XQT 
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Er- ,.,.  The  fy.tem  z^zzf^ssrz  sssu  SoSsSri; 


@ERS 

@ASG,T  8..F2 

@FOR,N  MEDICA*BIOCHEM. EXPFITPLOT  TPFS 

@MAP,N 

IN  TPF$. 

@XQT 

@ADD  MEDI CA*EXPFITPLOT . EXP2 854 
@COPY  8. ,MEDICA*INHIBPLOT. 

0MSG.W  PLOT , JRLOWE , X3836 
@PLOT , N MEDICA*INHIBPLOT. 

Plotting  progMm^^shSn^riec^ion  Hr"  dUrlng  ******  °peration  of  the 
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B.2.  Flowchart  for  Program  EXPFITPLOT. 
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B.3.  FORTRAN  Statements  for  Program  EXPFITPLOT 


C ME0ICA*3I0CHEH. EXPFITPLOT 


130  FORMAT  (/.0X.T1< 


/ I 


110  FORMAT  (2A6) 

123  FORMAT  (A3> 

1 33  FORMAT  (A3*  2X4) 

143  FORMAT  (E15.81 
153  FORMAT  (E13.4) 

1 60  FORMAT  ( I 3> 

1 73  FORMAT  ( F8. 3»F9.4> 

183  FORMAT  (3X.T1.*  PLOT  PARAMETERS  FOR  *.2A6) 

193  FORMAT  (3X.T1.*  ».A2,*  ZERO  =*.F5.2*#  INCHES 

SN.  ASLEN  :».F5.2.*  IN.*) 

2.33  FORMAT  ( 3 X . T1  , • 

i • 1 

210  format  ( 3X*  T1 » • The  PLOTTING  file  has  been  completed. m 


ASVMP  =**F5.2*»  I 


DIMENSION  KFIEL0(1003» .TIME(202»  *RATE(  202! 

8RITE  (6*130) 

CALL  PLOTS  (KFIELD.  1303,8  ) 

CALL  PLOT  ( 0. 0,-36. 3*-3) 

CALL  PLOT  (0.0, 2. 0.-3) 

CALL  SYMBOL  ( 1 . 0* 1 3. 3 , 3 . 1 4 ,26H JR LO BE  8I0MED  LAB  EXT  3836*0.0*26) 
220  RE  AO  (5,1 13,END=403)  IDRUN»ID2 
IF  ( lORUN.ro.  »NONE  *).  GO  TO  400 
RE  AO  (5.120»END  = 403>  TUNIT 
RE  AO  (5,1 13,END  = 40a)  VUNIT 
READ  (5,133  *END  = 400.)  OBSK  , OB  SK  A,  OB  SK  B 
READ  ( 5,1 33,EN0  = 400)  T IMO . T I M3  A . T I MOB 
READ  (5,143*EN0=400)  TIM3N 
READ  (5,1 33»END  = 40D)  V INF , VI NF A, VI NF B 
READ  (5*153, EN0=403)  VINFN 
RE  AO  ( 5,1 43 ,END=400)  HALF 
00  390  J-  1 . 4 
READ  (5«163*EN0=430>  N 
IF  (N)  400.430*230 

230  READ  (5,173.END=400)  ( TI ME ( I ) * RA TE ( I ) , 1= 1 * N) 

READ  (5,163. END=403)  M 
IF  ( (J-?*( J/2) l.EQ.O)  GO  TO  300 
240  IF  ( J.EO.  1)  60  TO  253 
PLOTIO=*B.» 

RATE(N*1) =3.3 

RATE!N*2)=VINFN/6.3 

TIME(N*1I =TIMON 

TIME(N»2>  =( S.»HALF-TIMON  J/9.5 

GO  TO  260 

250  CALL  SCALE  < TIME • 9 . 5* N*«l  I 
CALL  SCALE  ( R A TE * 6 . 0, N * *1  ) 
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u u u o 


«>L0TI0  = * A.» 

263  IF  ( TUNIT .NE. * SEC*  ) GO  TO  270 

CALL  AXIS  ( 3.0.3.0.15HTIME  IN  S£  CONDS* -1  5*  9 . 5 * 0.  0*  TI  HE  (N«l ) • 
*TIME(N*2) ) 

SC  TO  293 

270  IF  (TUNIT.NE.’HIN* ) GO  TO  280 

CALL  AXIS  (0.0*0. 0.15HTI ME  IN  MINUTES*  -15*9. 5*0. 0*  TIME  (N  ♦!  I* 
ATIK£(N*2> ) 

GO  TO  293 

283  CALL  AXIS  (0.0*0.0*1SH  TIME  IN  HOURS  • -15*9. 5*0. 0* TIME (N*l  I • 
STIME(N*2I I 

293  CALL  SVMBOL  ( 3.  3*  8 . 6*  0. 21  » 1 6HE  IP  ONCNTI  AL  PL0T»0.0»16) 

GC  TO  313 

333  CALL  SYMBOL  ( 6. C • 2 . S. 3 . 1 4 . SH V : V. 0.0*5) 

CALL  SYMBOL  ( 999 . • 2 . 46* C . 07  * 4H INF  *0.0.4) 

CALL  SVMBOL  ( 999 . • 2 . 5 * 3 . 1 4 *6HC 1 - E. 0.0*6) 

CALL  SYMBOL  (999. *2.6*0. 38*1  H-  «0.C*1  ) 

CALL  SYMBOL  (999.* 999. ,0.053*1 HK *0.0*1  ) 

CALL  «»LOT  (7. 9. 2. 6. 3) 

CALL  PLOT  (7.9*2.69 .?) 

CALL  PLCT  17.9*2.6*2) 

CALL  SYMBOL  ( 999 . » 2 . 58  * 3 . 04 » 3H 09 S. 0. 0* 3) 

CALL  SYMBOL  ( 999 . . 2 . 6 . 0 . 3 8 .4 H( T- T. 0. 0* 4 ) 

CALL  SYMBOL  ( 999 . * 2 . 5 8 . 3 . 04 « 1H0. 0. 0* 1) 

CALL  SYMBOL  ( 999 . * 2. 6 . 3 . 3 8 *1 H) *0 .0 . 1 ) 

CALL  SVMBOL  ( 999  . * 2 . 5 . 3 . 1 4 .1  H3  *0  .0  * 1 ) 

CALL  SYMBOL  ( 6. Q • 1 . 6 • 3 . 1 4 * 1H V • 0. 0* 1 ) 

CALL  SYMBOL  ( 999  . • 1 . 56  »Q  . 07  * 4H  INF  *0.0*4) 

CALL  SYMBOL  ( 999 . • 1 . 6 * 3. 1 4 * 1 H= *Q  .Q *1  ) 

CALL  SYMBOL  ( 999 . * 999. * 3 . 1 4 . VI NF *0 .0 . 3 ) 

CALL  SYMBOL  ( 999 . . 999.  . 0 . 1 4 . VI NF A . 0. 0* 4 ) 

CALL  SYMBOL  ( 999 . • 1 . 6 2.0 . 1 0 » 3H  X *0.0*3) 

CALL  SVMBOL  ( 999 . • 1 . 6 * 3 . 1 4 *2H1 0. 0. 0* 2) 

CALL  SYMBOL  ( 999 . . 1 . 7 . 3 . 3 8 . V INFB *0 .0 *4 ) 

CALL  SVMBOL  ( 6. 0 * 1 . 2* 0 .09  * 1 HK. 0. 0 * 1 ) 

CALL  PLOT  (6. 3*1. 2*3) 

CALL  PLOT  (6.0*1.35*2) 

CALL  PLOT  16.0*1.2*2) 

CALL  SYMBOL  ( 6. 09.  1 . 1 B .3 .07 . 4H  08  S *0.0*4) 

CALL  SYMBOL  ( 6. 42* 1 . 2 * C . 1 4 . 1 H= *0  .0  * 1 ) 

CALL  SYM90L  ( 999. * 999. *0. 14 , 09SK *0 .0  * 3 ) 

CALL  SYM30L  ( 999 . . 999. • 3 . 14  * 09 SK A , 0. 0. 4 ) 

CALL  SYMBOL  ( 999 .* 1. 2 2.0 . 1 3 . 3H  X ,0.0*3) 

CALL  SYMBOL  ( 999. » 1 . 2 . 3. 1 4 *2 HI  0. 0. 0. 2) 

CALL  SYMBOL  ( 999 . * 1 . 3 * 3 . 3 8 *OBS K3 *0  .0 *4  ) 

CALL  SVMBOL  ( 6. 0 *0 . 9* 0 . 14 , 1 HT , 0. 0 » 1 ) 

CALL  SYMBOL  ( 999 . * 0 . 78 . 0 . 07 . 4H0  *0.0*4) 

CALL  SYMBOL  ( 999 . * 0. 8 . 0. 1 4 * 1 Hr ,0  .0 » 1) 

CALL  SYMBOL  ( 999 . . 999. * 0. 1 4 . TI MO .0 .0 . 3 I 
CALL  SYM8CL  ( 999. » 999.  *C . 1 4 , TI MO A . 0. 0* 4 ) 

CALL  SYMBOL  ( 999 . * 3. 8 2 *0 . 13 . 3H  X *0.0*3) 
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CALL  SYMBOL  ( 999. • 3. 8 * 3 . 1 4 .2H1 3. 0. 0* 2) 

CALL  SYMBOL  I 999 . . 3. 9. 0 . 3 8 .T IHOB  .0  .0  . 4 ) 

TIME(N*l>  zT0 
TIM£(N*2)rTS 
RATE(N.1>=RB 
RAT£(N«2> =RS 

CALL  LINE  (TIME.RATE.N.l.M.U) 

GO  TO  320 

31G  CALL  L INC  (TIME.RAT£,N.l.M.O) 

CALL  SYMBOL  I 4 . 0 . 6 . 4 . 0 . I A * 8H  RUN-  10  I .0  .0  .8  ) 

CALL  SYMBOL  ( 5.  1 » 6 . 4 • 0 . 14  . IORUN.  0.  0.6  ) 

CALL  SYMBOL  ( 999 . . 394. , 0 . 1 4 . 102 * 0. 0. 6 I 
CALL  °L0 T (0.C»0.0«3) 

TB:TIMein*I ) 

TS=TIM£|n*2» 

3BZRAT£CN*1 ) 

rs=rate(n*2) 

50  TO  393 

323  2FRC  = -TIMEcn*  n/TIME  (N.2) 

CALL  PLCT  ( 7ER0.0.3.-3) 

IF  ( VUNIT .NE. • X RE  AC  * I GO  TO  33G 

CALL  AXIS  ( 3. 3.0.0. 14H*  RE AC TI VA TI ON . 1 4 t G . 0* 90 .0 .R AT E ( N* 1 ) 
iRATE<N«2>  > 

3C  TO  363 

333  IF  C VUNIT. NE. ’ABSORB*)  GC  TO  340 

CALL  AXIS  ( 0.0*0. 3> 13HA8S0R3ANCE  *1 0*6. 0.90»Q.RATE(  N* II • 
«RATE( N«2I ) 

GO  TO  363 

343  IF  ( VUN I T . NE . • x ACTI’)  GO  TO  350 

CALL  AXIS  (0.0.0.0.13H*  ACT  I VI  TV .1 0* 6 . 0. 90 .0 »R ATE ( N» 1 ) » 
iRAT£(N*2>>  , 

GO  TO  360 

353  CALL  AXIS  I 0. 0. 0 . 0 . 8HVEL0CIT Y * 0* 6 . 0.90 .0 »R ATE  I N«- 1 ) . 
iRAT£(N*2> ) 

363  ASYMPzVINFN/RATE (N*2  ) 

ASLEN=TB/TS»9.S 
IF  (ASVMP.GT.6.31)  GO  TO  370 
CALL  PLOT  < 0.0*ASYMP. 3 ) 

CALL  PLCT  ( ASLEN#ASYMP#2) 

370  CALL  PLOT  I 15 .3 . 0 . 3 . - 3 ) 

IF  (J.E0.4)  GO  TO  380 
WRITE  (6*180)  IDRUN* 102 
380  WRITE  (6.190)  PLOT  ID.  2ER0.  ASVMP.  ASLEN 
390  CONTINUE 

WRITE  (6.230) 

GO  TO  223 

403  CALL  PLOT  (15.0.3.3.949) 

WRITE  (6.210) 

WRITE  (6.130) 

STOP 

END 
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RUN-ID:  EXP2854 


8567 


RUN-ID:  EXP2854 


B 


•2 


-0.86  16.44  33.74  51.04  68.34  85.64  102.94  120.23  137.53  154.83 

TIME  IN  HOURS 


RUN-ID:  EXP120376 
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12.00  16.00  20.00  24.00  28.00  32.00 

TIME  IN  MINUTES 


RUN-ID:  EXP  1 20376 


TINE  IN  MINUTES 


B.5.  FORTRAN  Statements  for  Program  PLOTTEREXPFT 


Mf  DICA'SnCHEM 


PLOTTEREXPFT 


133 

FORMAT 

( / .3X.T 

113 

CCPKAT 

nx.  ti  • 

• ENTER, TOUR  NAME  USING  SIX 

characters 

123 

FORMAT 

( AS) 

1 3C 

FORMAT 

( 3 X • Tl  . 

• WHAT  IS  THE  OATA-IO?») 

1 4C 

FORMAT 

( 3X. Tl * 

•aERS* » 

1 50 

FORMAT 

13X.T1. 

•SASG. T 8 . .F2  * ) 

163 

FORMAT 

1 3X.  Tl  . 

•aFOR.N  MtDICA  *9  IOCHEM.E  XPFI  TPLOT.  TPF* 

1 73 

FORMAT 

( 3 X » T ! . 

•SNAP. N* ) 

1 8J 

FORMAT 

(3X.  Tl  . 

•IN  TP  Ft . • ) 

1 S3 

FORMAT 

( 3X»  Tl  . 

•L  13  M IS  0*  RL  0T  . • ) 

203 

FORMAT 

( OX. Tl  . 

• NEXT  OAT  A- ID?  * ) 

210 

F0PMAT 

( I A** ) 

223 

format 

( 3X. Tl  , 

•axciT* ) 

233 

FORMAT 

(3X.  TI  , 

•aAOO  MEDICAwEXPFI  TPICT.  ».  3A4) 

2 43 

FORMAT 

( 3X. Tl  , 

•aCOPY  8.  .MEOICA  wINHIBPLOT  .• 

) 

250 

FORMAT 

( OX. Tl  , 

•3MSG.W  PLOT ,* »A6. ».X3836* ) 

263 

FORMAT 

( 3X. Tl  , 

•aPLOT.N  MEDICA* INHI3PL0T. • ) 

270 

format 

( 3X. Tl  , 

283 

FORMA’- 

( 3 X.  Tl  , 

• NOW  TYPE  IN  «3  ADD  T£MO.«») 

WRITE 

(6.  130) 

WRITF 

(6.110) 

PFAO  (5.123)  NAME 

W 0 I T F 

(6. 130) 

wRITF 

(16.143) 

WRITE 

( 16.153) 

WRITE 

(16.163) 

WRITE 

( 16.1 73  ) 

WRITE 

( 16.183) 

WRITE 

( 15.1 93) 

30  TO 

303 

2 90 

WRITE 

( 6.  230 

300 

READ  (5.210)  I DR UN <102.103 

IF  (IDRUN.E0.* 

NONE’)  GO  TO  310 

WRITF 

( 16.223) 

W°  IT  E 

( 16.233) 

IDRUN.IC2.I03 

WRITE 

(16.243) 

WRITE 

< 16.253  ) 

NAME 

wD  I T E 

( 16.263) 

WRITE 

(6.270) 

30  TO 

293 

310 

WRITE 

(6.130) 

WPITE 

(6. 2E3) 

WRITE 

(6.130) 

WRITF 
5 TO® 
END 

( 16.143) 

OP  LESS.*) 
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B.6.  Sample  of  Teletype  Interactive  Output  for  Program  EXPFITPLOT 


#5ITE-ID 

*EDGEMOOD  ARSENAL  1108  SYSTEM  UER.  31-244  UPD  D(RSI)* 
SPUN  JRLOWE ? COST-CODE « PROJECT- 1 D > 30 - 1 000 
DATE:  G62176  TIME:  170818 

>@abd  Med  i cr*plotempf  i t . 

FURPUR  0026-06-2 1-17: 08 

READY 

RERDY 

FOR  S0E3-86- 21''76-17:08:46  (0«  ) 

END  FOR 

MRP28R1  RL71-3  06/21. -"76  17:08:50 

EMU  MRP 

ENTER  YOUR  MHME  USING  SIX  CHARACTERS  OR  LESS. 

: JRLOIIE 

UHRT  IS  THE  DATA- ID? 

>E::P8854 


Ht  XT  DATA- ID? 
ru  HIE 

MOM  TYRE  IN  "@RDD  TEMP." 

•'-:X:««3K%3KaSK»JKiK»«*^K3KIKK3KIK3KiKiKjK*jKm«»!iKSKi«K3KIK3l£%3K3RKBHKIK3IE)KKKKIKKiKKiKMm«!IHRilHR)KXiKK)nKMII( 
0RDD  TEMP. 

FURPUR  0026-06- '21 -17: 10 
READY 

F OP  S0E3-06-  2 1 -'76- 17:  10: 39  i' O « i 
END  FOR 

MAP20R1  RL7 1-3  86/21/76  17:10:50 
END  MAP 

,-\i:^j»*w>:};:a*;a^>:f»:«**r«*«*'femw*.*M***^«**i(i*»:*!*******i«**!i!i***** 
PLOT  PARAMETERS  FOR  FXP2854 

A.  FRO  •••  8.00  INCHES  ASYMP  = 5.00  IN.  ASLEM  - 7. SO  IN. 

X ZERO  « .05  INCHES  RSYMP  « 6.00  IN.  ASLEM  » 9. 45  IN. 

THE  PLOTTING  FILE  HRS  BEEN  COMPLETED. 

-se3e5*5«j(f*siofj«3*fsesj*5-!(tj*f5*eafae50S3OES*jSK3*f50f|O8^f5*E3af3oS5if5*t5*55lf*S»E**l»f3*EmS5»tmg»f*lieSiE»E3l6»63*E3iES»tl»tSlf3»e3K»eSlfllCIKI»6meaie*69it»HK«aitllt»t 

FI  If  PUR  0026-06/2 1-17:11 
i DLOCf  COPIED 
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DISTRIBUTION  LIST  2 


Names 

EDGEWOOD  ARSENAL 

CHIEF,  DEM1LTTARIZATION/DISPOSAL  OFFICE 
Attn:  SAREA-DM-E/Dr.  Ward 
CHIEF,  LEGAL  OFFICE 
CHIEF,  SAFETY  OFFICE 
CHIEF,  PLANS  OFFICE 
Attn:  SAREA-PL 
PUBLIC  HEALTH  SERVICE  LO 
AUTHOR'S  COPY,  Biomedical  Laboratory 

BIOMEDICAL  LABORATORY 
Attn:  SAREA-BL-M 
Attn:  SAREA-BL-R 
Attn:  SAREA-BL-RM 
Attn.  SAREA-BL-V 

CHEMICAL  LABORATORY 
Attn:  SAREA-CL 
Attn:  SAREA-CL-B 
Attn:  SAREA-CL-BS 
Attn:  SAREA-CL-D 
Attn:  SAREA-CL-P 
Attn:  SAREA-CL-TE 

DEVELOPMENT  AND  ENGINEERING  DIRECTORATE 
Attn:  SAREA-DE-D 
Attn:  SAREA-DE-DB 

DIRECTOR  OF  PRODUCT  ASSURANCE 
Attn:  SARF.A-PA-A 
Attn:  SAREA-PA-Q 

DIRECTOR  OF  TECHNICAL  SUPPORT 
Attn:  SAREA-TS-R 
Attn:  SAREA-TS-L 

DEPARTMENT  OF  DEFENSE 

Administrator 

Defense  Documentation  Center 
Attn:  Accessions  Division 
Cameron  Station 
Alexandria,  VA  22314 

Director 

Defense  Intelligence  Agency 
Attn:  DB-4G1 
Washington,  DC  20301 

DEPARTMENT  OF  THE  ARMY 

HQDA  (DAMO-SSC) 

WASH  DC  20310 

Director 

Defense  Civil  Preparedness  Agency 
Attn:  RE 

Washington,  DC  20301 


Copies  Names 

Chief,  Office  of  Research, 

Development  & Acquisition 
Attn:  DAMA-ARZ-D 
1 Washington,  DC  20310 

1 

1 CINCUSAREUR 

Attn:  AEAMD-PM 
1 APO  New  York  09403 

8 

56  US  ARMY  HEALTH  SERVICE  COMMAND 

Commander 

1 US  Army  Environmental  Hygiene  Agency 

1 Attn:  HSE-AA-L,  Librarian,  Bldg  2100 

1 APG-Edgewood  Area 

1 

Superintendent 
Academy  of  Health  Sciences 
1 US  Army 

1 Attn:  HSA-CDC 

1 Attn:  HSA-IHE 

1 Fort  Sam  Houston,  TX  78234 

1 

1 US  ARMY  MATERIEL  DEVELOPMENT  AND 

READINESS  COMMAND 

1 Commander 

1 US  Army  Materiel  Development  and 

Readiness  Command 
Attn:  DRCLDC 
1 5001  Eisenhower  Ave 

1 Alexandria,  VA  22333 

Commander 

2 US  Army  Foreign  Science  & Technology  Center 

3 Attn:  DRXST-IS1 
220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

Commander 
DARCOM,  STITEUR 
12  Attn:  DRXST-STf 

Box  48,  APO  New  York  09710 

Commander 

US  Army  Science  & Technology 
Center-Far  East  Office 
1 APO  San  Francisco  96328 

US  ARMY  ARMAMENT  COMMAND 

Commander 

1 US  Army  Armament  Command 

Attn:  DRSAR-ASN 
Attn:  DRSAR-IMB-C 
Rock  Island,  IL  61201 

1 


Copies 


1 


1 


1 


1 

1 


1 


2 


1 


1 


1 

1 
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DISTRIBUTION  LIST  2 (Contd) 


Names  Copies 

Commander 

US  Army  Dugway  Proving  Ground 

Attn'  Technical  Library,  Docu  Sect  1 

Dugway,  UT  84022 

Commander 

Rocky  Mountain  Arsenal 

Attn:  SARRM-EA  1 

Denver,  CO  80240 

US  ARMY  TRAINING  & DOCTRINE  COMMAND 


Commandant 

US  Army  Missile  & Munitions 
Center  & School 

Attn:  ATSK-DT-MU-EOD  1 

Redstone  Arsenal,  AL  35809 

Commandant 

US  Army  Military  Police  School 

Attn:  ATZN-CD  1 

Fort  McClellan,  AL  36205 

Commander 

US  Army  Infantry  Center 

Attn:  ATZB-CD-MS-C  1 

Fort  Bcnning,  GA  31905 

Commandant 

US  Army  Ordnance  Center  & School 

Attn:  ATSL-CD-MS  1 

Attn:  ATSL-DTG  1 

APG-Abcrdeen  Area 


US  ARMY  TEST  & EVALUATION  COMMAND 

Record  Copy 
CDR,  APG 

Attn:  STEAP-AD-R/RHA  1 

APG-Edgewood  Area,  Bldg  E5179 

CDR, APG 

Attn:  STEAP-TL  1 

APG-Aberdeen  Area 

Commander 

US  Army  Tropic  Test  Center 

Attn:  STETC-MO-A  (Tech  Library)  1 

APO  New  York  09827 

Commander 

Dugway  Proving  Ground 

Attn:  STEDP-PO  1 

Dugway,  UT  84022 

DEPARTMENT  OF  THE  NAVY 


Names  Copies 

Commander 

Naval  Explosive  Ordnance  Disposal  Facility 

Attn:  Army  Chemical  Office  1 

Indian  Head,  MD  20640 

Commander  1 

Naval  Surface  Weapons  Center 

Attn:  Tech  Lib  & Info  Svcs  Br  1 

White  Oak  Laboratory 
Silver  Spring,  MD  20910 


Commander 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 

Attn:  DX-21  1 

Dahlgren,  VA  22448 

Commander 
Naval  Weapons  Center 

Attn:  Technical  Library/Code  233  1 

China  Lake,  CA  93555 

DEPARTMENT  OF  THE  AIR  FORCE 

HQ  Foreign  Technology  Division  (AFSO 

Attn:  PDRR  1 

Wright-Patterson  AFB.OH  45433 

Commander 

Aeronautical  Systems  Division 

Attn:  ASD/AELD  1 

Wright-Patterson  AFB.OH  45433 

Director 

Air  Force  Inspection  and  Safety  Center 

Attn:  IGD(AFISC/SEV)  1 

Norton  AFB.CA  92409 

Commander 

Armament  Development  & Test  Center 

Attn:  DLOSL  (Technical  Library)  1 

Eglin  AFB,  FL  32542 

OUTSIDE  AGENCIES 

Director  of  Toxicology  1 

National  Research  Council 
2101  Constitution  Ave,  NW 
Washington,  DC  20418 


Chief  of  Naval  Research 

Attn:  Code  443  1 

800  N.  Quincy  Street 
Arlington,  VA  22217 
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